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Abstract 

In the present paper we overview our recent results on intrinsic frictional prop- 
erties of adsorbed monolayers, composed of mobile hard-core particles undergoing 
continuous exchanges with a vapor phase. Within the framework of a dynamical 
master equation approach, describing the time evolution of the system, we determine 
in the most general form the terminal velocity of some biased impure molecule - the 
tracer particle (TP), constrained to move inside the adsorbed monolayer probing its 
frictional properties, define the frictional forces as well as the particles density distri- 
bution in the monolayer. Results for one-dimensional solid substrates, appropriate 
to adsorbtion on polymer chains, are compared against the Monte Carlo simulation 
data, which confirms our analytical predictions. 



1 Introduction. 



In the present paper we overview our recent results on intrinsic friction of monolayers 
emerging on solids exposed to a vapor phase |1]-§J. Such layers are involved in various 
technological and material processing operations, including, for instance, coating, gluing 
or lubrication. Knowledge of their intrinsic frictional properties is important for conceptual 
understanding of different transport processes taking place within molecular films, film's 
stability, as well as spreading of ultrathin liquid films on solid surfaces |||[], spontaneous 
or forced dewetting of monolayers @-10| or island formation on solid surfaces [11 . 



Since the early works of Langmuir, much effort has been invested in the analysis of the 
equilibrium properties of the adsorbed films |Il2|-|i6[| . Here, significant analytical results 
have been obtained predicting different phase transitions and ordering phenomena, which 
well agree with the available experimental data. As well, some approximate results have 
been obtained for both dynamics of an isolated adatom on a corrugated surface and collec- 
tive diffusion, describing spreading of the macroscopic density fluctuations in interacting 
adsorbates being in contact with the vapor |T7|-Pm • 

Another important aspect of dynamical behavior concerns tracer diffusion in adsorbates, 
which is observed experimentally in STM or field ion measurements and provides a useful 
information about adsorbate's viscosity or intrinsic friction. This problem is not only a 
challenging question in its own right due to emerging non-trivial, essentially cooperative 
behavior, but is also crucial for understanding of various dynamical processes taking place 
on solid surfaces. However, most of available theoretical studies of tracer diffusion in 



adsorbed layers (see, e.g. Refs.pT|-p9|) do exclude the possibility of particles exchanges 
with the vapor. 

Here we focus on this important issue and provide a theoretical description of the 
properties of tracer diffusion in adsorbed monolayers in contact with a vapor phase - a 
reservoir of particles. More specifically, the system we consider consists of (a) a solid 
substrate, which is modeled in a usual fashion as a regular lattice of adsorbtion sites; (b) a 
monolayer of adsorbed, mobile hard-core particles in contact with a vapor and (c) a single 
hard-core tracer particle (TP). We suppose that the monolayer particles move randomly 
along the lattice by performing symmetric hopping motion between the neighboring lattice 
sites, which process is constrained by mutual hard-core interactions, and may desorb from 
and adsorb onto the lattice from the vapor with some prescribed rates dependent on the 
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vapor pressure, temperature and the interactions with the solid substrate. In contrast, the 
tracer particle is constrained to move along the lattice only, (i.e. it can not desorb to the 
vapor), and is subject to a constant external force of an arbitrary magnitude E. Hence, 
the TP performs a biased random walk, constrained by the hard-core interactions with 
the monolayer particles, and always remains within the monolayer, probing its frictional 
properties. 

The questions we address here are the following: First, focussing on one- and two- 
dimensional systems, we aim to determine the force- velocity relation. That is, as t — ► oo, 
the TP ultimately attains a constant velocity V tr (E), which depends on the magnitude of 
the applied external force; the functional form of this dependence in the most general case 
constitutes the primary goal of our analysis. Next, we study the form of the force- velocity 
relation in the limit of a vanishingly small external bias. This allows us, in particular, 
to show that the frictional force exerted on the TP by the monolayer particles is viscous, 
and to evaluate the corresponding friction coefficient. Lastly, we analyse how the biased 
TP perturbs the particles density distribution in the monolayer. As a matter of fact, we 
proceed to show that there are stationary density profiles around the TP as t — > oo, which 
mirror a remarkable cooperative behavior. 

We finally remark that our model can be viewed from a somewhat different perspective. 
Namely, on the one hand, the model under study is a certain generalization of the "tracer 
diffusion in a hard-core lattice gas" problem (see, e.g. Ref.[ZT] for an extensive review) to 
the case where the random walk performed by the TP is biased and the number of particles 
in the monolayer is not explicitly conserved, due to exchanges with the reservoir. We recall 
that even this, by now classic model constitutes a many-body problem for which no exact 
general expression of the tracer diffusion coefficient D tr is known. On the other hand, our 
model represents a novel example of the so called "dynamical percolation" models |5H| - 



35], invoked to describe transport processes in many situations with dynamical disorder. 
In this context, the particles of the monolayer can be thought of as representing some 
fluctuating environment, which hinders the motion of an impure molecule - the TP, which 
might be, for example, a charge carrier. An important aspect of our model, which makes 
it different from the previously proposed models of dynamic percolative environments, is 
that we include the hard-core interaction between "environment" particles and the tracer 
particle, such that the latter may influence itself the dynamics of the environment. Lastly, 
we note that the model under study can be thought of as some simplified picture of the 
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stagnant layers emerging in liquids being in contact with a solid body. It is well known (see, 
e.g. Ref.||36||) that liquids in close vicinity of a solid interfaces - at distances within a few 
molecular diameters, do possess completely different physical properties compared to these 
of the bulk phase. In this "stagnant" region, in which an intrinsically disordered liquid 
phase is spanned by and contends with the ordering potential of the solid, liquid's viscosity 
is drastically enhanced and transport processes, (regarding, say, biased diffusion of charged 
carriers in solutions), are essentially hindered. Thus our model can be viewed as a two-level 
approximate model of this challenging physical system, in which the reservoir mimics the 
bulk fluid phase with very rapid transport, while the adsorbed monolayer represents the 
stagnant layer emerging on the solid-liquid interface. 

The paper is structured as follows: In Section |2| we formulate our model in case of a 
two-dimensional solid substrate and introduce basic notations. We write down then the 
dynamical equations which govern the time evolution of the monolayer particles and of the 
tracer, and outline the decoupling approximation used to close this system of equations. 
These equations are presented for the general, <i-dimensional case. Sections [I] and [| are 
respectively devoted to the presentation of the results of this general approach for one- 
dimensional substrates, which situation is appropriate to adsorbtion on polymer chains 
|37|, 38], and also for two-dimensional solid substrates. Finally, we conclude in Section |5| 
with a brief summary and discussion of our results. 



2 The model and basic equations 
2.1 The model. 

Consider a two-dimensional solid surface with some concentration of adsorbtion sites, which 
is brought in contact with a reservoir containing identic, electrically neutral particles - a 
vapor phase (Fig.|TJ) , maintained at a constant pressure. For simplicity of exposition, we 
assume here that adsorbtion sites form a regular square lattice of spacing a. We suppose 
next that the reservoir particles may adsorb onto any vacant adsorbtion site at a fixed rate 
f/r*, which rate depends on the vapor pressure and the energy gain due to the adsorbtion 
event. Further on, the adsorbed particles may move randomly along the lattice by hopping 
at a rate l/4r* to any of 4 neighboring adsorbtion sites, which process is constrained by 
hard-core exclusion preventing multiple occupancy of any of the sites. Lastly, the adsorbed 
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particles may desorb from the lattice back to the reservoir at rate g/r*, which is dependent 
on the barrier against desorption. Both / and g are site and environment independent. 




Figure 1: Adsorbed monolayer in contact with a vapor. Grey spheres denote the monolayer 
(vapor) particles; the smaller black sphere stands for the driven tracer particle. 

Note that such a model represents, of course, a certain idealization. In "real world" 
systems, the adsorbed particles experience two types of interactions: the solid-particle (SP) 
interactions, i.e. interactions with the atoms of the host solid, and the particle-particle 
(PP) interactions, i.e. mutual interactions with each other. Now, both the SP and the 
PP interactions are characterized by a harsh repulsion at short scales, and are attractive 
at longer distances. The SP repulsion keeps the adsorbed particles some distance apart 
of the solid, while the SP attraction favors adsorbtion and hinders particles desorption 
as well as migration along the solid surface. In this regard, our model corresponds to 
the regime of the so-called intermediate localized adsorption |I5|,[16|: the particles forming 
an adsorbed monolayer are neither completely fixed in the potential wells created by the 
SP interactions, nor completely mobile. This means, that the potential wells are rather 
deep with respect to particles desorption (desorption barrier ^> kT), such that the 
adsorbate forms a submonolayer, but nonetheless, have a much lower energy barrier Vi 
against the lateral movement across the surface, Ud 3> Vi > kT. In this regime, each 
adsorbed particle spends a considerable part of its time at the bottom of a potential well 
and jumps sometimes, solely due to the thermal activation, from one potential minimum 
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to another in its neighborhood; after the jump is performed, the particle dissipates all its 
energy to the host solid. Thus, on a macroscopic time scale the particles do not possess 
any velocity. The time r* separating two successive jump events, is just the typical time 
a given particle spends in a given well vibrating around its minimum; t* is related to the 
temperature, the barrier for lateral motion and the frequency of the solid atoms' vibrations 
by the Arrhenius formula. 

We emphasize that such a type of random motion is essentially different of the standard 
hydrodynamic picture of particles random motion in the two-dimensional "bulk" liquid 
phase, e.g. in free-standing liquid films, in which case there is a velocity distribution and 
spatially random motion results from the PP scattering. In this case, the dynamics may 
be only approximately considered as an activated hopping of particles, confined to some 
effective cells by the potential field of their neighbors, along a lattice-like structure of 
such cells (see, e.g. Refs. [j39| [40|1). In contrast to the dynamical model to be studied here, 
standard two-dimensional hydrodynamics pressumes that the particles do not interact with 
the underlying solid. In realistic systems, of course, both the particle-particle scattering 
and scattering by the potential wells due to the interactions with the host solid, (as well as 



the corresponding dissipation channels), are important [^q, [fljl. 111 particular, it has been 
shown that addition of a dissipation to the host solid removes the infrared divergencies in 
the dynamic density correlation functions and thus makes the transport coefficients finite 
P^.p3|. On the other hand, homogeneous adsorbed monolayers may only exist in systems 
in which the attractive part of the PP interaction potential is essentially weaker than that 
describing interactions with the solid; otherwise, such monolayers become unstable and 
dewet spontatneously from the solid surface. As a matter of fact, for stable homogeneous 
monolayers, the PP interactions are at least ten times weaker that the interactions with 
the solid atoms |nj . 

Consequently, the standard hydrodynamic picture of particles dynamics is inappropri- 
ate under the defined above physical conditions. Contrary to that, any adsorbed particle 
moves due to random hopping events, activated by chaotic vibrations of the solid atoms, 
along the local minima of an array of potential wells, created due to the interactions with 
the solid |15|,|16[|. As we have already remarked, in the physical conditions under which 
such a dynamics takes place, the PP interactions are much weaker than the SP interactions 
and hence do not perturb significantly the regular array of potential wells due to the SP 
interactions. In our model, we discard completely the attractive part of the PP interac- 
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tion potential and take into account only the repulsive one, which is approximated by an 
abrupt, hard-core-type potential. 

To describe the occupancy of lattice sites, we introduce a time-dependent variable i](TV), 
which may assume two values: 



r,(R) 



1, if the site R is occupied by an adsorbed particle, 
0, if the site R is empty. 



Note that the local variable r/(R) can change its value due to adsorption, desorption and 
random hopping events. Note also that random hopping events do conserve the total 
number of adsorbed particles and hence, the average adsorbate density p s (t). On the other 
hands, adsorption/desorption processes do change r](R) locally such that the total number 
of particles in the adsorbed monolayer is not explicitly conserved. However, the mean 
density of the adsorbate, p s (t) =< r?(R) >, approaches as t — > oo a constant value 

p ' = 7+9 (1) 

This relation is well known and represents the customary Langmuir adsorption isotherm 



12|| . We finally remark that in the analysis of the stationary-state behavior, we can always 
turn to the conserved particles number limit by setting / and g equal to zero and keeping 
their ratio fixed, i.e. supposing that f/g = p s /(l — p s )- This limit will correspond to the 
model of biased tracer diffusion in a hard-core lattice gas with fixed particles density p s , 
and will allow us to check our analytical predictions against some already known results 



21^p5|, p9| . We will furnish such a comparison at the end of the section f|. 

Further on, at t = we introduce at the lattice origin an extra hard-core particle, whose 
motion we would like to follow; position of this particle at time t is denoted as R tr . Note 
that the tracer particle (TP) is designed to measure the resistance offered by the monolayer 
particles to the external perturbance, or, in other words, to measure the intrinsic frictional 
properties of the adsorbate. 

Now, we stipulate that the TP is different from the adsorbed particles in two aspects: 
first, it can not desorb from the lattice and second, it is subject to some external driving 
force, which favors its jumps into a preferential direction. Physically, such a situation may 
be realized, for instance, if this only particle is charged and the system is subject to a 
uniform electric field E. We suppose here, for simplicity of exposition, that the external 
force E is oriented according to the unit vector ei. 
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The dynamics of the biased TP is defined in the usual fashion: We suppose that the TP, 
which occupies the site H tr at time t, waits an exponentially distributed time with mean 
r, and then attempts to hop onto one of 4 neighboring sites, Tt tr + e^, where e„ are 4 unit 
vectors of the hypercubic lattice. In what follows we adopt the notation v = {±1,±2}, 
where ±ex will denote the direction of the external force E. Next, the jump direction is 
chosen according to the probablity p u , which obeys: 



where (3 is the reciprocal temperature, (E ■ e) stands for the scalar product, the charge of 
the TP is set equal to unity and the sum with the subscript \i denotes summation over all 
possible orientations of the vector e M ; that is, /i = {±1, ±2}. 

After the jump direction is chosen, the TP attempts to hop onto the target site. The 
hop is instantaneously fulfilled if the target site is vacant at this moment of time; otherwise, 
i.e., if the target site is occupied by any adsorbed particle, the jump is rejected and the 
TP remains at its position. 

2.2 Evolution equations. 

Now, we derive the evolution equations in a general, d- dimensional case, which will allow 
us to compare the behavior emerging in one- and two-dimensional systems. We begin by 
introducing some auxiliary definitions. Let 77 = {r/(R)} denote the entire set of the occu- 
pation variables, which defines the instantaneous configuration of the adsorbed particles 
at the lattice at time moment t. Next, let P(R tr , ?7;0 stand for the joint probability of 
finding at time t the TP at the site R tr and all adsorbed particles in the configuration 77. 
Then, denoting as rf' v a configuration obtained from 77 by the Kawasaki-type exchange of 
the occupation variables of two neighboring sites r and r + e„, and as if - a configuration 
obtained from the original 77 by the replacement r/(r) — > 1 — 77(1*), which corresponds to 
the Glauber-type flip of the occupation variable due to the adsorption/desorption events, 
we have that the time evolution of the configuration probability P(R tr , 77; t) obeys the 



exp | 




(2) 



Pu = 



f(E-e M r 
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following master equation: 

d 

d t P(R tr , V ;t) = —J2 E { P ( Rtr ' ^ ') " P ( Rtr ' ^ £ ) } 

A*=l r^R tr -e M ,R tr 

+ - E^{ (1 " ^(Rtr)) P(Rtr - e M , 77; t) - (1 - ?7(Rtr + 0) P ( R tr, W *) } 

+ ^E {(l-^(r))P(Rt r ^ r ;t)-^(r)P(R tr ,r/;t)} 

r^Rtr 

+ ^E {^(r) J P(Rt r ,^;t)-(l-77(r))P(R tr ,r ? ;t)}. (3) 

The mean velocity Vt r (t) of the TP can be obtained by multiplying both sides of Eq.(|3|) 
by (Rtr • ei) and summing over all possible configurations (Rt r , 77). This results in the 
following exact equation determining the TP velocity: 

V »® = J t S (Rtr ' ei ) P ( Rtr ' K t) = ^{pi (l - Hex; *)) - p_! (l - fc(e_ i; £)) }, (4) 

Rtr,?? 

where 

fc(A; t) = E ^( R tr + A)P(Rtr, ^; *) (5) 

Rtr,?) 

is the probability of having at time t an adsorbed particle at position A, defined in the 
frame of reference moving with the TP. In other words, fc(A; t) can be thought of as being 
the density profile in the adsorbed monolayer as seen from the moving TP. 

Equation @ signifies that the velocity of the TP is dependent on the monolayer par- 
ticles density in the immediate vicinity of the tracer. If the monolayer is perfectly stirred, 
or, in other words, if k(X; t) = p s everywhere, (which implies immediate decoupling of Rt r 
and 77), one would obtain from Eq.@ a trivial mean-field result 

^ 0) = (Pi-P-i)(1-Ps)-, (6) 

r 

which states that the only effect of the medium on the TP dynamics is that its jump time 
r is merely renormalized by a factor (1 — p s )~ 1 , which represents the inverse concentration 
of voids in the monolayer; note that then (1 — p s )/r defines simply the mean frequency of 
successful jump events. 

However, the situation appears to be more complicated and, as we proceed to show, 
k(X;t) is different from the equilibrium value p s everywhere, except for |A| — > 00. This 
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means that the TP strongly perturbs the particles distribution in the monolayer - it is no 
longer uniform and some non-trivial stationary density profiles emerge. 

Now, in order to calculate the instantaneous mean velocity of the TP we have to 
determine the mean particles density at the neighboring to the TP sites R tr + e±i, which 
requires, in turn, computation of the density profile k(X; t) for arbitrary A. The latter can 
be found from the master equation (|3|) by multiplying both sides by 7/ (Rtr) an d performing 
the summation over all configurations (R tr ,r/). In doing so, we find that these equations 
are not closed with respect to k(\;t), but are coupled to the third-order correlations, 

T(X, e„, t)=J2 ^( R tr + AMR*.- + e /J p (Rtr, V) *) (?) 

Rtr,»7 

In turn, if we proceed further to the third-order correlations, we find that these are coupled 
respectively to the fourth-order correlations. Consequently, in order to compute V tr , one 
faces the problem of solving an infinite hierarchy of coupled equations for the correlation 
functions. Here we resort to the simplest non-trivial closure of the hierarchy in terms of 
k(X;t), which has been first proposed in Refs.[26] and [27], and represent T(A, e u ;t) as 



vC&tr + A)r/(R tr + e At )P(R tr , rj; t) 



Rtr,»7 



« ( J2 ^ R *r + A ) P ( R tr, V] t) ) ( J2 ^( Rtr + e M) P (^tr, V, t) 

= k(X;t)k(e^t), (8) 

Some arguments justifying such an approximation a posteriori are presented in Sections 3 
and 4 (see also Ref. 

Using the approximation in Eq.(|8]), we obtain 

2dr*d t k(X; t) = Lk(X; t) + 2df, (9) 

which holds for all A, except for A = {0, ±e 1; e 2 . . . , e d }. One the other hand, for these 
special sites A = e u with v = {±1, 2, . . . , d} we find 

2dr*d t k{e u - t) = {L + A u )k(e„- 1) + 4/, (10) 

where L is the operator 

L^^A M V M -4(/ + ^), (11) 
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and the coefficients are defined by 

2dr* 

A^t) = 1 + — P M (1 - fc(e M ; t)). (12) 

Note that Eq.([T0|) represents, from the mathematical point of view, the boundary condi- 
tions for the general evolution equation (^), imposed on the sites in the immediate vicinity 
of the TP. Equations (0) and QUID together with Eq. (|) thus consitute a closed system of 
equations which suffice computation of all properties of interest. 

2.3 Stationary solution of the evolution equations 

We turn to the limit t — > oo and suppose that both the density profiles and stationary 
velocity of the TP have non-trivial stationary values 

k(X) = lim fc(A;t), V tr = lim V tr (t), and A, ee lim AJt) (13) 

Define next the local deviations of k(\) from the unperturbed density as 

h{\) EE k(X) - Ps (14) 

Choosing that h(0) = 0, we obtain then the following fundamental system of equations: 

Lh(X) = for all A ^ {0, e±i, . . . , e ±d }, (15) 

{L + A u )h(e„) + Ps (A u - A_ u ) = for A = {0, e ±1 , . . . , e ±d }, (16) 

which determine the deviation from the unperturbed density p s in the stationary state. 
Note also that in virtue of an evident symmetry, h(e u ) = h(e_ u ) and A v = A_ v for 

ue{2,... ,d}. 

The general approach to solution of coupled non-linear Eqs.(|),(|l5]) and ([16|) has been 
discussed in detail in Ref. [§J. Here we merely note that despite the fact that using the 
decoupling scheme in Eq.(^) we effectively close the system of equations on the level of 
the pair correlations, solution of Eqs.QTB]) and ( |T6"D (or, equivalently, of Eqs.(^) and (]TTJ|) ) 



still poses serious technical difficulties: Namely, these equations are non-linear with respect 
to the TP velocity, which enters the gradient term on the rhs of the evolution equations 
for the pair correlation, and does depend itself on the values of the monolayer particles 
densities in the immediate vicinity of the TP. Solution of this system of non-linear equation 
for one- and two-dimensional substrates is displayed in two next sections. 



10 



3 One-dimensional adsorbed monolayer 



For one- dimensional lattices, which which situation is appropriate to adsorbtion on polymer 
chains |37|,[3£|, general solution of Eqs. fll5|) and ( fUf) has the following form: 



k n = k(X) = p s + K±exp(^ — <r\n\/\±\ , A = an, n E Z, 



(17) 



where the characteristic lengths A± obey 



A- 



=F a In 1 



A l + A_! + 2(f + g) t J (A x + A_! + 2(/ + <?)) - 4A^_ 



while the amplitudes K± are given respectively by 



A_i — Ai exp(— cr/A H 



and 



(19) 



(20) 



'A_iexp(-cr/A_) - Ai 

Note that A_ > A+, and consequently, the local density past the TP approaches its non- 
perturbed value p s slower than in front of it; this signifies that correlations between the TP 
position and particle distribution are stronger past the TP. Next, K + is always positive, 
while K_ < 0; this means that the density profile is a non-monotoneous function of A and 
is characterized by a jammed region in front of the TP, in which the local density is higher 
than p s , and a depleted region past the TP in which the density is lower than p s . 

Now, we are in position to obtain a system of two closed-form non-linear equations 
determing implicitly the unknown parameters Ai and A_i, which will allow us to compute 
the TP terminal velocity, related to A±i through V tr = cr(Ai — A_i)/2r*. Substituting 
Eq.flTFD into Eq.(0), we find 

Ai - A_! 



Ai = 1 + 



and 



A_ 



1 + 



Pit' 



T 



P-i-r 



T 



1 - p s - p s 



A_iexp(a/A' + ) - A x 



(21) 



Ps - Ps 



Ai - A_i 



A_! — A 1 exp(a/A / _ 



(22) 
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Figure 2: Terminal velocity of the probe molecule as a function of the adsorption probabil- 
ity / at different values of the parameter g. The probe hopping probabilities are pi = 0.6 
and p-i = 0.4. The solid lines give the analytical solution while the filled squares de- 
note the results of Monte-Carlo simulations. Upper curves correspond to g — 0.8, the 
intermediate - to g — 0.5 and the lower - to g — 0.3, respectively. 

Resolution of this system leads then to the stationary velocity of the TP (see Fig.|]) as well 
as the density profiles (see Fig.|3|). For arbitrary values of p, f and g the parameters A±±, 
defined by Eqs. (pT|) and (0), and consequently, the terminal velocity V tr can be determined 
only numerically (see Figs. 2 to 4). However, V tr can be found analytically in the explicit 
form in the limit of a vanishingly small force E, E — > 0. Expanding A±i in the Taylor 
series in powers of E and retaining only linear with E terms, we find that the TP velocity 
follows 

V tr ~ C'E, (23) 

which relation can be thought off as the analog of the Stokes formula for driven motion 
in a one- dimensional adsorbed monolayer undergoing continuous particles exchanges with 
the vapor phase. Equation (p3|) signifies that the frictional force exerted on the TP by the 
monolayer particles is viscous. The friction coefficient, i.e. the proportionality factor in 
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Eq.(p3"D is given explicitly by 



2r 



C = — i -i ?±L z (9 1) 

f3a 2 (l-p s )l r(f + g) 1 + ^1 + 2(1 + r*(l - p s )M/(f + gY 

Note that the friction coefficient in Eq. (|24|) can be written down as the sum of two contri- 
butions ( = ( cm + (coop- The first one, ( cm = 2r/ !3a 2 {l — p s ) is a typicial mean-field result 
and corresponds to a perfertly homogeneous monolayer (see discussion following Eq.([5|)). 
The second one, 

8r*p s 1 
Cc ° OP " ^ 2 (1 - Ps)(f + 9) 1 + y/1 + 2(1 + r*(l - p a )/r)/(/ + ^ (25) 

has, however, a more complicated origin. Namely, it reflects a cooperative behavior emerg- 
ing in the monolayer, associated with the formation of inhomogeneous density profiles (see 
FigfJ) - the formation of a "traffic jam" in front of the TP and a "depleted" region past the 
TP (for more details, see The characteristic lengths of these two regions as well as the 
amplitudes K± depend on the magnitude of the TP velocity; on the other hand, the TP 
velocity is itself dependent on the density profiles, in virtue of Eq.(£|). This results in an 
intricate interplay between the jamming effect of the TP and smoothening of the created 
inhomogeneities by diffusive processes. Note also that cooperative behavior becomes most 
prominent in the conserved particle number limit P^ ETTR . Setting f,g—>0, while keeping 
their ratio fixed (which insures that p s stays constant), one notices that ( COO p gets infinitely 



large. As a matter of fact, as it has been shown in Refs. |26| and in such a situation 
no stationary density profiles around the TP exist; the size of both the "traffic jam" and 



depleted regions grow in proportion to the TP mean displacement X tr (t) ~ \ft [26 



Consequently, in the conserved particle number limit ( COO p grows indefinitely in proportion 
to Vi- 
lli order to check our analytical predictions, we have performed numerical Monte Carlo 
simulations of the exact Master equation of the problem, using the method of Gillespie 
p7| . Results of these simulations, performed at different values of the parameters /, g, 
and pi, are also represented in Figs.^] and |^. 

Consider finally the situation with E = 0, in which case the terminal velocity van- 
ishes and one for which one expects conventional diffusive motion with the mean square 
displacement of the form 



Xlit) = 2D tr t, (26) 
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Figure 3: Density profile around stationary moving probe molecule for / = 0.1, g = 0.3 
and p = 0.98. The solid line is the plot of the analytical solution. Filled squares are the 
results of Monte-Carlo simulations. 



where D tr is some unknown function of the system parameters. Heuristically, we can 
compute D tr for the system under study if we assume the validity of the Einstein relation 
Dt r — ft/C between the friction coefficient and the self-diffusion coefficient D tr of the TP 
48fl , which yields 



D 



tr 



a\l - Ps 
2t 



1 



<f + 9)1 + ^l + 2(l + T*(l-p a )/T)/(f + g) 



(27) 



tr 



Monte Carlo simulations (see Fig.[|) of the system evidently confirm our prediction for D 
given by Eq.(p7|), and hence, confirm the validity of the Einstein relation for the system 
under study. This is, of course, not an unexpected, but still a non-trivial result 



4 Two-dimensional adsorbed monolayer 

We turn now to the case of a two-dimensional substrate - adsorbtion onto the surface of a 
solid exposed to a vapor phase. Here, the situation gets somewhat more difficult from the 
computational point of view; we have now to solve the partial difference equations problem 
(|T5D, ( P^D rather than the mere difference equations arising in the one-dimensional case. 
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0.0 0.5 1.0 

/ 

Figure 4: Self-diffusion coefficient of the probe molecule as a function of the adsorption 
probability /. Notations and values of g are the same as in Figs.0. 

Solution in two-dimensions can be found in a most convenient fashion if we introduce 
the generating function for the particle density profiles, defined as 

+oo +oo 

H{w u w 2 )= E ^,n 2 «, (28) 

ni=— CO U2 =— CO 

where h nijTl2 = h{niex + n 2 e 2 ). Multiplying both sides of Eqs. (|i~5|) and (^) by w'^w^ 2 , 
and performing summations over n\ and n 2 , we find that H(wi,W2) is given explicitly by 

H(w 1 ,w 2 ) = -K(wx,w 2 ) {AiWi 1 + A^wi + A 2 {w 2 + w^ 1 ) - a} , (29) 

where a = A u + 4(/ + g) and 

K(w u w 2 ) = J2 Mw\u\ Ul ~ l)He u ) + PsiAt - A^)( Wl - wf 1 ). (30) 

Equations (|29| ) and (|30| ) determine the generation function for the density profiles exactly. 

Before we proceed to the inversion of H(wi, w 2 ) with respect to the variables w± and w 2 , 
we note that we can access interesting integral characteristic of the density profiles directly 
using the result in Eqs. (|29|) and (|30D . Namely, as we have already remarked, the presence 
of the driven TP induces an inhomogeneous density distribution in the monolayer. One 
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can thus pose a natural question whether equilibrium between adsorption and desorption 
processes gets shifted due to such a perturbancy, i.e. whether the equilibrium density in the 
monolayer is different from that given by Eq. ([!]). The answer is trivially "no" in the case 
when the particles number is explicitly conserved, but in the general case with arbitrary 
/ and g this is not at all evident: similarly to the behavior in one-dimensional system 
one expects that also in two-dimensions the density profiles are asymmetric as seen from 
the stationary moving TP and are characterized by a condensed, "traffic-jam" -like region 
in front of and a depleted region past the TP. One anticipates then that the desorption 
events are favored in front of the TP, while the adsorption events are evidently suppressed 
by the excess density. On the other hand, past the TP desorption is diminished due to the 
particles depletion while adsorption may proceed more readily due to the same reason. It is 
thus not at all clear a priori whether these two effects can compensate each other exactly, 
in view of a possible asymmetry of the density profiles, as it happens in the one-dimensional 
model (see Fig.|3|). 

For this purpose, we study the behavior of the integral deviation Q of the density 
from the equilibrium value p s , i.e. f2 = ^n^-oo Sn°^-oo ^m,n 2 ) which can be computed 
straightforwardly from Eqs. (|29|) and (0) by setting both w\ and w 2 equal to unity. Noticing 
that K{wi — 1,W2 = 1) = 0, and that A\ + A_i + 2A 2 — a = — 4(/ + g), i.e. is strictly 
negative as soon as adsorption/desorption processes are present, we obtain then that Q is 
stricly equal to 0. This implies, in turn, that the perturbancy of the density distribution 
in the monolayer created by the driven TP does not shift the global balance between the 
adsorption and desorption events. 

Inversion of the generating function with respect to w± and w 2 requires quite an involved 
mathematical analysis, which has been presented in detail in Ref. M. General solution for 
the density profiles reads: 

K u n 2 = a- l ^2A v h{e u )V- v F num - p s (A 1 - A^)(Vt - V_i)F nii „ 2 }. (31) 

V 

with 

/ A ,\ ni / 2 r°° / \ 

F num = l-^J e-*I m (2a- 1 y/A^tj I n2 (2a~ 1 A 2 t) dt, (32) 

where l n (z) stands for the modified Bessel function. We mention that F ni n2 has an inter- 
esting physical interpretation in terms of the generating function of a random walk of a 
single particle (that is, in absence of the particles environment) ||44|| . 
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Now, the Eqs.^D and (0) display h m,n 2 as & function of the coefficients A v that remain 
to be determined. As a matter of fact, these coefficients depend themselves on the local 
densities in the immediate vicinity of the tracer, i.e. on h(e u ). This implies that we have 
to determine them from Eqs.(|3lD and ( |32"D in a self-consistent way ||. Some analysis (see 
0] shows that A v are determined implicitly as the solution of the following system of three 
non-linear matrix equations 



At* 

{±1,2}, A u = l + p p 

T 



1 - p s - p s (A 1 - A-i 



detq 
' det (7 



(33) 



where 



C 



(A 1 V- 1 F ei -a A_xViF ei A 2 V- 2 F ei \ 

AiV-iF^ A^ 1 V l F e _ 1 - a A 2 V. 2 F e _ 1 
y AiV_iF ea A_xViF e2 A 2 V„ 2 F e2 -aJ 



(34) 



the matrix C v stands for the matrix obtained from C by replacing the z/-th column by the 
column-vector F, 



1 (Vi - V_x)F ei ^ 
F = (V! - V-i)^ , 

\ (V! - V_!)F e2 , 

while the local deviations h(e v ) are expressed in terms of A v as 

%„) = ( 



(35) 



Ps) + 



T 



1 - A, 



Lastly, the TP terminal velocity obeys 



~(pi 1 + p 



4r* px det Ci - p_i det C_i 



-i 



(36) 



(37) 



r - i r det C 

which represents the desired general force- velocity relation for the system under study, 
which is valid for arbitrary magnitude of the external bias and arbitrary values of other 
system's parameters. 



4.1 Asymptotical behavior of the density profiles at large sepa- 
rations from the TP. 

Asymptotical behavior of the density profiles at large distances from the TP follows from 



the analysis of the analyticity properties of the complex function N(z) = J2 n= _ 00 h n flZ n . 
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It has been shown in Ref. M that in front of the TP, the deviation h n ^ always decays 
exponentially with the distance: 



exp ^ — n/\j^j 



h n,0 ~ K + -1/2 » ( 38 ) 



where the characteristic length A + obeys: 

A + = In" 1 f-^ <( ^ - AH J(^-A 2 ) -A 1 A_ 1 }\ (30) 




Note that A + stays finite for any values of the system parameters. 

On contrary, the behavior of the density profiles at large distances past the tracer 
qualitatively depends on the physical situation studied. In the general case when exchanges 
with the particles reservoir are taken into allowed, the decay of the density profiles is still 
exponential with the distance: 



exp ^ — n/X^j 



h. nfl ~ K_ ^— (40) 

where 

A - = - ln-1 (i{f-W(^) 2 -^-]) (41) 

Note that in the general case the characteristic lenghts again, similarly to the one- dimensional 
case, satisfy the inequality A_ > A + , which means that the correlations between the TP 
and the particles of the monolayer are always stronger past than in front of the TP. 

Such correlations can even become extremely strong in the special case when the par- 
ticles exhanges with the vapor phase are forbidden, i.e. in the conserved particles number 
limit, which can be realized for the monolayers sandwiched in a narrow gap between two 
solid surfaces. In this case, we have that A_ becomes infinitely large and, in the limit 
n — > +oo, the deviation of the particle density from the equilibrium value p s follows 

K'_ ( 3 / 1 

h- n ,0 = m (l + — +0[ — 

Remarkably enough, in this case the correlations between the TP position and the particles 
distribution vanish algebraically slow with the distance! This implies, in turn, that in the 
conserved particles number case, the mixing of the monolayer is not efficient enough to 
prevent the appearence of the quasi-long-range order and the medium "remembers" the 
passage of the TP on a long time and space scale, which signifies very strong memory 
effects. 



l + Z+ b - ( 4 2) 
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4.2 Limit of small applied force 

We turn now to the limit f3E -C 1, in which case the problem simplifies considerably and 
allows to obtain explicit results for the local densities in the immediate vicinity of the TP 
and consequently, for the TP terminal velocity and its diffusivity. 

In this limit, we arrive again at a Stokes-type formula of the form Vt r ~ E/Q-, where 
now 

(3a\l -Ps) y r (f + g + 1 + r * (1 _ ps y T \ _ x ' 



with 



and 



2f + g+l + T*(l-p s )/T> l44J 



oo 

-t 



C(x) = U e-^(I (ari) - h(xt))I (xt)dt j . (45) 

Note that we again are able to sigle out two physical meaningful contributions to the friction 
coefficient (. Namely, the first term on the rhs of Eq.(^3|) is just the mean-field-type result 
corresponding to a perfectly stirred monolayer, in which correlations between the TP and 
the monolayer particles are discarded. The second term, similarly to the one-dimensional 
case, mirrors the cooperative behavior emerging in the monolayer and is associated with 
the backflow effects. In contrast to the one-dimensional case, however, the contribution to 
the overall friction coefficient stemming out of the cooperative effects remains finite in the 
conserved particles limit. 

We also wish to remark, that a qualitatively similar physical effect has been predicted 
recently for a different model system involving a charged particle moving at a constant 
speed a small distance above the surface of an incompressible, infinitely deep liquid. It 
has been shown in Refs. |f45|, [46fl , that the interactions between the moving particle and the 



fluid molecules induce an effective frictional force exerted on the particle, producing a local 
distortion of the liquid interface, - a bump, which travels together with the particle and 
increases effectively its mass. The mass of the bump, which is analogous to the jammed 
region appearing in our model, depends itself on the particle's velocity resulting in a non- 
linear coupling between the medium-induced frictional force exerted on the particle and 
its velocity [JK], ^ . 
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Lastly, assuming a priori that the Einstein relation holds for the system under study, 
we estimate the TP diffusion coefficient D tr = f3~ 1 (~ 1 as 

Ar = — (1 - Pa) S 1 



4r v riy [ r 4(/ + ^+l + r*(l-p.)/r)£(x)-l + (3p a -l)r*/r 

(46) 

It seems now interesting to compare our general result in Eq.(fi6|) against the classical 
result of Nakazato and Kitahara ||22]| , which describes TP diffusion coefficient in a two- 
dimensional lattice gas with conserved particles number. Setting / and g equal to zero, 
while assuming that their ratio has a fixed value, / / g = p s j (I — p s ), we have then that 

D tr = ^(1 -Ps) {l - PS T 4(1 - r , (1 _ Ps)/r)C (l/2) - 1 + (3p s - l)r* /r } • (47) 



£(1/2) = , - r , (48) 



Using next the equality |49|: 



we find that the right-hand-side of Eq. (|47|) attains the form 



a 2 



D ><- - T^ 1 ~ Pi I 1 " — TTn w n irvTTTi . w S > ^ ( l9) 

4r [ T 1 + I 1 - Ps)t*/t - (1 - 2/7r)(l + (1 - 3p s )T*/r) 



which expression coincides exactly with the earlier result obtained in Refs.[^] and |23 
within the framework of a different, compared to ours, analytical techniques. The result 
in Eq.(|45|) is known to be exact in the limits p s < 1 and p s ~ 1, and serves as a very 
good approximation for the self-diffusion coefficient in hard-core lattice gases of arbitrary 
density [pi]], which supports in a way the validity of the approximation invoked in Eq.(P). 



5 Conclusion 

To conclude, we have studied analytically the intrinsic frictional properties of adsorbed 
monolayers, composed of mobile hard-core particles undergoing continuous exchanges with 
the vapor. Our analytical approach has been based on the master equation, describing 
the time evolution of the system, which has allowed us to evaluate a system of coupled 
dynamical equations for the TP velocity and a hierarchy of correlation functions. To solve 
these coupled equations, we have invoked an approximate closure scheme based on the de- 
composition of the third-order correlation functions into a product of pairwise correlations, 
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which has been first introduced in Ref.f26|| for a related model of a driven tracer particle 
dynamics in a one- dimensional lattice gas with conserved particles number. Within the 
framework of this approximation, we have derived a system of coupled, discrete-space equa- 
tions describing evolution of the density profiles in the adsorbed monolayer, as seen from 
the moving tracer, and its velocity V tr . We have shown that the density profile around the 
tracer is strongly inhomogeneous: the local density of the adsorbed particles in front of 
the tracer is higher than the average and approaches the average value as an exponential 
function of the distance from the tracer. On the other hand, past the tracer the local den- 
sity is always lower than the average, and depending on whether the number of particles is 
explicitly conserved or not, the local density past the tracer may tend to the average value 
either as an exponential or even as an algebraic function of the distance. The latter reveals 
especially strong memory effects and strong correlations between the particle distribution 
in the environment and the carrier position. Next, we have derived a general force- velocity 
relation, which defines the terminal velocity of the tracer particle for arbitrary applied 
fields and arbitrary values of other system parameters. We have demonstrated next that 
in the limit of a vanishingly small external bias this relation attains a simple, but physi- 
cally meaningful form of the Stokes formula, which signifies that in this limit the frictional 
force exerted on the tracer by the adsorbed monolayer particles is viscous. Correspond- 
ing friction coefficient has been also explicitly determined. In addition, we estimated the 
self- diffusion coefficient of the tracer in the absence of the field and showed that it reduces 



to the well-know result of Refs. [P^1 and |23[ in the limit when the particles number is 
conserved. 
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